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Abstract 

The slow dynamics of a system as it approaches a phase transition, associated 
with the slowing down in the decay of a correlation function, can be caused by 
a sharp increase in the probability of a particle's returning to its original state 
following a transition, rather than to a slowing down in the transition rates as 
is usually assumed. The results of our calculations show that this is the case 
for the ferromagnetic Potts model. The implications of this result for various 
theories of the glass transition are discussed. 
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I. Introduction 



The drastic slowing down of a system's dynamics as the temperature T 
is lowered towards a critical temperature T c is a well-known feature of many 
physical systems. Typical examples include the slowing down of the a-relaxation 
process in supercooled liquids as the glass transition temperature is approached 
and the critical slowing down of the dynamical processes in a system undergoing 
a phase transition as the phase transition temperature is approached. There are 
two aspects of this slowing down, namely the general theory and the nature of 
the specific mechanism. For the general theory, in the case of phase transitions 
this slowing down can readily be derived quite simply from the Landau theory 
of phase transitions, which involves an expression for the free energy of the 
system in terms of an order parameter and its fluctuations. P while for the glass 
transition there are various suggestions for general theories but none that is yet 
generally acceptedPl With regard to the nature of the specific mechanism, it 
is usually assumed that this slowing down is associated with a drastic slowing 
down in the rates at which the particles make transitions. In this paper, we 
show that there is another possible mechanism, namely a change in the results 
of successive transitions, and present results for a specific model system where 
this is the case. 

The calculation of a system's properties is usually performed in terms of 
correlation functions C(t), that describe the average relationship between the 
states of the particles at time t with their original state at time t = Q. One reason 
for this is that, according to Kubo's formula, many of a system's observable 
macroscopic properties are determined by correlation functions. For instance, 
the frequency dependent dielectric susceptibility x( w ) is proportional to the 
Fourier transform of the dipole-dipole correlation function < M(t).M(0) >P1 
Hence, a rapid increase in the probability that molecules return to their original 
state (in the relevant phase space) following a transition, which we call the 
return probability, will also lead to a drastic slowing down in the rate of decay 
of the correlation functions, and so of the associated relaxation processes and 
of the system's dynamics. The question of whether this process, rather than a 
drastic slowing down of the transition rates, is responsible for the slow dynamics 
is well worth investigating. For instance, for supercooled liquids it has important 
implications for the type of theory that can explain the glass transition, as we 
discuss below. Incidentally, this type of approach is in line with that of Berthier 
and GarrahanP who emphasized the importance of considering transitions in 
real space. 

The property of a system that is most often studied, both theoretically and 
experimentally, is the correlation function, but this on its own cannot show 
which of the above mechanisms is responsible for the slowing down. Another 
property that can readily be studied theoretically (although not so easily mea- 
sured experimentally), and which can throw some light on this question, is the 
fraction of particles in the system that have never changed their state up to time 
t. We call this function the system's transition function and denote it by Pt(t). 
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If the particles move randomly after making a transition, the correlation and 
transition functions should decay with time in a very similar manner. However, 
if the particles have a tendency to return to their original state because of the 
arrangement of their neighbors, as in Funke's mismatch and relaxation model 
for ionic conductor^ and in the cage effect in mode coupling theory^ for in- 
stance, the correlation function will decay much more slowly with time than the 
transition function. It follows that a comparison of the change in the time de- 
pendence of these functions as the critical temperature is approached can throw 
light on which of the above mechanisms is responsible for the slow dynamics. 
In this paper, we make such a comparison for a simple model system in which 
there is a phase transition, and so critical slowing down as it is approached, and 
then consider the implications of our results for other systems. 

The system that we consider is the ferromagnetic g-state Potts model, for 
which we have recently reported results for the transition function Pt(t), in a 
papeiP henceforth referred to as I. The function Pt(t) was there called P(t) and 
termed the relaxation function, but we have changed its name and symbol in 
this paper so as to avoid any impression that it is related to a physical property 
such as the dielectric polarization. We have extended the calculations in I 
by calculating simultaneously the transition function Pt(t) and the normalized 
correlation function Cr(t), which we define below, as T c /T increases towards 
unity. In section 2, we describe briefly the model system and how we analyzed 
its properties. The results of our calculations are presented and discussed in 
section 3, while our conclusions are summarized in section 4. 

II. The Potts model and its analysis. 

The Hamiltonian for the ordered ferromagnetic q-spin Potts model with in- 
teractions only between the spins at adjacent sites can be written aJ^ 

where J > 0, the first sum is over all the sites i in the system and the second 
one over all the sites that are nearest neighbors of the site i, the spins 
(Tj can take any integer value between 1 and q, and 5 is the Kronecker delta, 
5(a, b) = 1 if a = b, if a ^ b. For this system, some care is needed in defining 
the relevant correlation function. The elementary correlation function is just 
s(t) =< 5[<7i(t), ai(0)] >, but this will never tend to zero as t — > oo since 
the spins can only assume q distinct values, so that even for completely random 
<Ji(t) a fraction l/q of them will equal <7i(0). Accordingly for a system of N 
spins we define the normalized correlation function Cr(t) by 

Cr{t) =< £>Mi), <*(())] - l/q} > /[N(l - l/q)] (2) 

i 

for which Cr(0) = 1 and if eventually there is no correlation between the initial 
and final states then Cr(oo) = 0. 
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It follows from equation (1) that the energy E associated with a site on 
which the spin is equal to that on z adjacent sites is just — z J. The probability 
of a change in the spin at a site which involves an increase of AE in this energy 
at temperature T was taken to have the standard form 

w = w , AE < (3) 
w = w exp{-AE/k B T), AE > 

For an infinite square lattice, the values of the critical temperature T c are given 
by the solutions v — v c of the equation v 2 — q, where v = exp(J/[fcsT]) — 1, 
while the phase transition is of second order for q < 4 and of first order for 
q > 4P In order to examine both types of phase transition, we performed our 
calculations for two systems, with values of q — 3, 6. For our simulations we 
considered, as in I, a square lattice of 200 x 200 sites, and we here chose Wo = 1 
since we are interested in temperatures close to T c . The simulation techniques 
used were the same as in I, apart from starting the anneals from an initial state 
of all identical spins instead of one with random spins. Once a steady state was 
reached, five sets of several successive simulation runs were performed on the 
system, with each run proceeding until the spin had changed at least once at 
99.9% of the sites. A maximum of 5 runs was performed in each set, but these 
were stopped once Cr{t) became less than 0.002. Finally, the average values of 
Pt(t) and Cr(t) for the five sets of runs were calculated. 

The results were analyzed in terms of the mean relaxation times, < r >, for 
the decay of these functions, and two different techniques were used to calculate 
them. In general, a monotonically decreasing function f(t) can be expressed 
as a sum or integral of functions cxp(— t/r), with a distribution g(r) of relax- 
ation times r, f(t) = J g(r) exp(— f/r)dr. In that case, < r >= J g(r)TdT = 
Jo f(t)dt, and this latter expression was used to calculate < r >. The main 
method that we used for our calculations was to fit Pt(t) and Cr(t), for val- 
ues of these functions greater than 0.01, to stretched exponential functions 
f(t) = Aexp[-(t/r ) /3 ], for which f™ f(t)dt = (A/P)T(l/j3)r . In all cases, 
an excellent fit was obtained (in a graph of ln[/(t)] as a function of t) with 
A = 1. We restricted our calculations to values of x = T c /T such that after a 
set of 5 runs Cr(t) < 0.1, since for higher final values of Cr(t) it is difficult to 
justify a fit to such a function. In addition, < t > was calculated directly for 
the two functions from J °° f(t)dt, and the value obtained was found to be very 
close to that obtained from the fit to the stretched exponential function (except 
for Cr(t) at the highest values of x, where the cut-off in Cr(t) after five runs led 
to a much lower value for J f(t)dt since the integral did not extend to negligible 
values of Cr(t) ). Accordingly, all the results for < r > that we present, which 
are denoted by Tp for Pt{t) and by tq for Cr(t), are those derived from fits to 
a stretched exponential function. 
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III. Results and Discussion. 



In figure 1, we show both ln(rp) and ln(rc) as functions of T c /T for q = 3 
and for q — 6. Incidentally, the results presented here for ln(rp) differ slightly 
from those in our previous paper I because there we showed the stretched ex- 
ponential relaxation time r rather than < t >, and also used w = 0.5. It is 
immediately apparent that while for T c /T — 0.85 the relaxation times tq and 
Tp are similar, the relaxation times tq for the correlation function becomes 
larger than those for the transition function Tp as T c /T increases, until for 
T c /T = 0.99 we find that tc/tp = 33 for q = 3 and 6.3 for q = 6. In order 
to show the difference between the functions Cr(t) and Pt(t), and not just the 
ratio of their relaxation times, we present these functions for T c /T — 0.85 and 
T c /T — 0.99 for q — 3 in figure 2 and for q = 6 in figure 3. These figures show 
clearly that for T c /T = 0.85 the transition and correlation functions decay with 
time in a very similar manner, but for T c /T = 0.99 the correlation function 
decays much more slowly than the transition function. We also note from figure 
1 that as T c /T — > 1 while ln(rp) increases fairly steadily, with no particular sign 
of a phase transition, ln(r c ) starts to increase very rapidly, so that for higher 
values of T c /T than those shown in the figure it decays very little on the time 
scale of around 5rp to which we extended our simulations. Thus, it is the rapid 
increase in tq (rather than in Tp) which is a sign of the critical slowing down 
associated with a phase transition in our system, and this is due to the rapid 
increase in the return probability. 

For the Potts model, there is a simple qualitative explanation of this increase 
in the return probability as the temperature decreases, and also of the difference 
between the systems with q = 3 and those with q = 6. For a site within a cluster 
of identical spins on a square lattice, the activation energy AE required for a 
transition is 4J. As a result, the slowing down of the decay of the transition 
function, which is associated with the temperature dependence of the factor 
exp(— AE/kpT) in equation (3) and with the increasing fraction of sites within 
clusters (as found in I) as T decreases, does not change dramatically as T 
approches T c . On the other hand, when the spin at a site within a cluster 
changes, it can change again without any activation energy until it returns to 
its original state, when it again requires an activation energy of 4J to change it. 
Thus, the spins on sites within clusters have a strong tendency to return to their 
original values unless the spins on the adjacent states have changed before they 
do so, something that happens most easily at higher temperatures, and also 
for small clusters since the spins at sites on the boundary of a cluster require 
a lower activation energy to change. As a result, the increase in the fraction 
of sites within clusters (and in the mean cluster size) as T decreases leads to 
a rapid increase in the return probability, and so in the ratio of tc/tp, with 
a dramatic increase as the transition to a single phase is approached. This is 
just an example of critical slowing down as the phase transition is approached. 
It sets in at lower temperatures (higher values of T c /T) for q = 6 than for 
q = 3 because the fraction of sites within clusters at a given value of T c /T 
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decreases with increasing number of different possible values of the spins, i.e. 
with increasing q, as was found in figure 8 of I. 

The above results show that for the ferromagnetic Potts model the main 
mechanism responsible for the critical slowing down as the phase transition is 
approached is the increasing return probability and not a drastic slowing down 
in the rate at which transitions occur. Such a mechanism is not restricted to 
the Potts model, and not necessarily even to phase transitions, and can well be 
relevant to the glass transition in supercooled liquids. As we discussed in I, the 
Potts model can be regarded as a simple model for plastic crystals, which have 
relaxation properties similar to those of supercooled liquidsPSl Accordingly, we 
now consider the relevance of this type of process to the drastic slowing down of 
the relaxation rates in supercooled liquids as the glass transition is approached. 
Our results show that this need not be caused by a change in the basic relaxation 
mechanism, and in fact no such change was found in recent experiments by 
Huang and RichertPH Instead, as for the Potts model, it can be due to the local 
molecular arrangements being such that even when a molecule instantaneously 
changes its position or orientation it rapidly returns to its original lower energy 
state. Such a process is obviously compatible with theories of the glass transition 
based on heterogeneities^ and/or smeared out phase transitions^! since the 
Potts model exhibits a phase transition while the clusters of states of identical 
spin in it are local heterogeneities. It is also fully compatible with theories of 
the glass transition that involve potential energy landscapes,^ for which the 
molecular transitions correspond to jumps within a meta-basin and those which 
lead to diffusion of the molecule (and so to decay of the correlation function) to 
jumps of the system from one meta-basin to another. However, such a process 
does not fit so well to theories based on free volume^or on kinetic constraints j^l 
where it is the basic transition rates of the molecules that change drastically 
as the temperature is lowered. Of course, it is quite possible that different 
mechanisms are the dominant ones in different systems. 

IV. Conclusions 

A comparison between the transition function Pt{t) and the correlation func- 
tion Cr{t) of a system that exhibits a drastic slowing down of its dynamical 
properties can show whether this slow dynamics is caused mainly by a slowing 
down in the transition rates of the particles or by an increase in their tendency 
to return to their original state following such a transition. For the ferromag- 
netic Potts model, the latter process was found to be responsible for the critical 
slowing down as the phase transition was approached. The question of which is 
the dominant process has important consequences for the different theories of 
the glass transition in various types of supercooled liquids. 
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Captions for figures 

Figure 1. The mean relaxation time r for the relaxation function Pt(t) (open 
symbols) and for the correlation function Cr(t) (closed symbols) of the g-state 
Potts model as functions of T c /T. The results for q = 3 are denoted by black 
squares, and those for q = 6 by red circles, while the broken lines are just to 
guide the eye. 

Figure 2. The relaxation function Pt(t) (continuous black lines) and the 
correlation function Cr{t) (broken red lines) for the 3-state Potts model, for 
T c /T = 0.85 (lower pair) and for T c /T = 0.99 (upper pair). 

Figure 3. The relaxation function Pt(t) (continuous black lines) and the 
correlation function Cr{t) (broken red lines) for the 6-state Potts model, for 
T c /T = 0.85 (lower pair) and for T c /T = 0.99 (upper pair). 
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